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Noise induced currents and reliability of transport in frictional ratchets 
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We study the coherence of transport of an overdamped Brownian particle in frictional ratchet 
system in the presence of external Gaussian white noise fluctuations. The analytical expressions 
for the particle velocity and diffusion coefficient are derived for this system and the reliability or 
coherence of transport is analysed by means of their ratio in terms of a dimensionless Peclet number. 
We show that the coherence in the transport can be enhanced or degraded depending sensitively on 
the frictional profile with respect to the underlying potential. 
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I. INTRODUCTION 

The study of the interplay of noise and nonlinear dy- 
namics in systems under nonequilibrium conditions has 
generated wide interdisciplinary interests in the last 
two decades. Noise or fluctuations, which are normally 
considered to be an hinderance, are found to play an 
active constructive role in nonequilibrium systems. In- 
tact, this constructive role of noise (as opposed to the 
conventional wisdom of its destructive or its disorgan- 
ising role) has become a new paradigm in the study 
of complex systems. In the so called ratchet systems 
the presence of spatial/temporal anisotropy in poten- 
tial together with nonequilibrium perturbations enable 
the extraction of useful work from random fluctuations 
without the violation of the second law of thermody- 
namics 0, 0- In such systems it is possible to in- 
duce directed motion from nonequilibrium fluctuations 
in the absence of bias. Much of the studies in differ- 
ent classes of ratchet models deal with the nature of 
currents and their reversals IllL stochastic energetics 
(thermodynamic efficiency) |3, 1^ etc. However, trans- 
port of Brownian particle is always accompanied by a 
diffusive spread and this spread is intimately related to 
the question of reliability or quality of transport. The 
diffusive spread infact detriments the quality of trans- 
port. There exists very few studies which address the 
question of diffusion accompanying transport in ratchet 
systems 0, IE ■ o^'' present work we address this 
aspect of transport and study the coherence in trans- 
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port in a frictional ratchet in the presence of an exter- 
nal parametric Gaussian white noise fluctuation. This 
is infact studied in terms of a dimensionless quantity 
called the Peclet number (Pe) which is the ratio of ve- 
locity to the diffusion constant. Higher the Pe, lesser is 
the diffusive spread and higher is the transport coher- 
ence. Infact, subcellular transport in biological systems 
amidst a noisy environment are modeled based on the 
principle of ratchet mechanism and the experimental 
studies on these molecular motors show them to have 
highly efficient and reliable transport with Peclet num- 
ber ranging from 2 to 6 . A value of Pe greater than 
2 corresponds to coherent transport 0. The Peclet 
numbers for some of the models like flashing and rock- 
ing ratchets were found to be ~ 0.2 and ~ 0.6 fi^ re- 
spectively implying a less reliable transport. Another 
study on symmetric periodic potentials along with spa- 
tially modulated white noise showed a coherent trans- 
port with Peclet number less than 3. In the same study 
a special kind of strongly asymmetric potential is found 
to increase Pe to 20 in some range of physical param- 
eters Idi. 



There exists many physical systems like flashing 
ratchets 0, rocking ratchets 0|, time asymmetric 
ratchets etc., where different aspects of noise induced 
transport has been widely studied 1, 2]. In the above 
models, to generate unidirectional current the nonequi- 
librium fluctuations need to be correlated in time. 
There exists a possibility to get unidirectional cur- 
rent even in presence of symmetric ratchet potentials 
provided it is driven by a time correlated asymmetric 
force In our present work we consider yet an- 

other class of ratchets, namely the frictional ratchets, 
where the friction coefficient and subsequently the dif- 
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fusion coefficient varies in space 0, 0| . In these fric- 
tional ratchets it is possible to get unidirectional cur- 
rents even in a symmetric underlying potential but in 
presence of an external noise which need not be corre- 
lated in time unlike in earlier models. 

Space dependent diffusion coefficient, D{x), felt by 
the Brownian particle could arise either due to space 
dependent temperature or space dependent friction co- 
efficient. In the frictional ratchet system which we 
consider in the present work the unidirectional current 
arises due to a combination of both space dependent 
friction coefficient and external parametric white noise. 
The temperature of the bath (or the environment) of 
the Brownian particle is characterised by a constant 
temperature T. In the presence of external paramet- 
ric noise the overdamped Brownian particle on an av- 
erage absorbs energy from the external noise source. 
The strength of the absorbed energy depends on the 
local frictional coefficient. Hence the problem of par- 
ticle motion in an inhomogeneous medium in presence 
of an external noise becomes equivalent to t he p rob- 
lem in a space dependent temperature 0, llsL . 
Such systems are known to generate unidirectional 
currents. This follows as a corollary to Landauer's 
blow torch theorem that the notion of stability changes 
latically 
ities Iq. 

valid in equilibrium systems, does not hold. 

Frictional inhomogeneities are common in supcrlat- 
tice structures, semiconductors or motion in porous 
media. Particles moving close to a surface experience 
space dependent friction j^l . It is believed that molec- 
ular motor protiens moving close along the periodic 
structures of microtubules experience a space depen- 
dent friction (17i |. Frictional inhomogeneity changes 
the dynamics of the particle nontrivially as compared 
to the homogeneous case. This in turn has been shown 
to give rise to many counter intuitive phenomena like 
noise induced stability, stochastic resonance, enhance- 
ment in efficiency etc., in driven non-equilibrium sys- 
tems IHEl- 

In our present work we show that system inhomo- 
geneities may help in enhancing/degrading the coher- 
ence in the transport depending sensitively on the 
physical parameters. We emphasize mainly the case 
where the underlying potential is a simple sinusoidal 
symmetric potential. The role of spatial asymmetry in 
potential is also discussed. The external noise is found 
to play a constructive role in enhancing the coherence 



dramatically in the presence of temperature inhomo 
geneities 15]. In such cases the notion of local stability. 



in transport. As opposed to this, temperature (internal 
fluctuations) degrades the coherence in transport. 



II. MODEL: 

We start with the Kramer's equation of motion for a 
Brownian particle of unit mass in contact with a heat 
bath in a medium with spatially varying friction coef- 
ficient ri{q) at temperature T. In addition an external 
parametric Gaussian white noise fluctuation ^ (t) is also 
included. The equation of motion is given by 



q = -v{q)q-v'{q) + ./k^T^f{t) + ^{t) (1) 

where V{q) is the potential seen by the Brownian par- 
ticle and f(t) is an internal Gausian white noise fluc- 
tuation arising from the bath having the property that 

< f{t) >^ 0, and < f{t)f{t') 2S{t - t') where 

< ... > denotes the ensemble average and q the co- 
ordinate of the particle. Also, < £,{t) >~ and 

< ^{t)^it') >= 2r(5(t - t'), where F is the strength 
of the external white noise C(^)- The above equation, 
Eq. ^ has been derived earlier from the microscopic 
consideration of system bath coupling 12, IS] . 

On time scales larger than the inverse friction coef- 
ficient, one can in most practical cases consider 
the overdamped limit of the Langevin equation. This 
in turn correspond to the adiabatic elimination of the 
fast variable, velocity, from the equation of motion by 
putting p = q = for a homogeneous system. In con- 
trast, for the case of inhomogeneous system the above 
method of elimination does not work and Sancho et 
al. lo'l has given a proper prescription for the elimina- 
tion of fast variables. The corresponding overdamped 
Langevin equation for the Brownian particle in a space 
dependent frictional medium is given by 



q = 



V'iq) kBTi{q) ^ ksT 
Viq) 2[rj{q)f V Viq) 



fit) 



viq)' 



(2) 



Using van Kampen Lemma ^20j and Novikov's the- 
orem |21| we get the corresponding Fokker-Planck 
or Smoluchowski equation for the probabilit y d ensity 
P{q, t) of a particle being at g at a time t as (2^ 



dP 
'dt 



d_ 

dq 



V'{q) U{q) 

i]{q) 773(g) 



P 



viq) 



r 



dP' 
rf{q)\ dq 
^ (3) 

For periodic functions V{q) and "qiq) with periodic- 
ity L a finite probability current is obtained and one 
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readily gets the analytical expression for the particle 
velocity as 



v = L- 



(l-exp(-5)) 



(4) 



Here '4>{q) is the dimensionless generalized effective po- 
tential given by 

j r][x)[kBT-q(x) +T\ 

and A{q) is the effective space dependent diffusion co- 
efficient given by 



A{q) 



kBTrjjq) + r 
7^2 (q) 



Also, 



S = ijj{q) -'4>{q + 2'K) 



(6) 



(7) 



determines the effective slope of the generalized poten- 
tial i^iq). Thus the sign of 5 gives the direction of 
current in Eq. ^ 

In our present work we have taken the potential 
V{q) = Vosin(g) and ■q{q) = ?7o[l - Asin(g - 0)], 
< A < 1. The phase lag </> between V{q) and ■q{q) 
brings in the intrinsic asymmetry in the dynamics of 
the system. 

For the case with T = we get a simple analytical 
expression for the effective potential as 



i>{q) = 



Vb?7o / A2:sin(0) A 



r V 2 

-I- sin(2a;) sin(</))] -f sin(a;) 



-[cos(2x)cos(0) (8) 



with 5 given by 



6 = - 



yo?/o7i'Asin((/)) 



ln[r7(x)] (9) 



(10) 



The first term in the right hand side of Eq. |^ 
represents the tilt in the effective potential. This 
tilt identically vanishes when </> = or tt. Hence 
it is expected that the unidirectional current does 
not arise for the case when the phase lag is or tt. 
From Eq. ^ and Eq. ^| it is clear that for < < tt 
current will be in the negative direction while for 
TT < < 27r the current will be in the positive 
direction. Fig. ^ shows the plot of the effective 
potential as a function of coordinate q for a fixed 
noise strength F = 0.36 for two different values of 




FIG. 1: Effective potential '(/'(g) as a function of q for Vb = 
1, r = 0.36 and A = 0.9 at zero temperature for (j) = O.SQvr 
and 0.897r. 



phase lag cj) between the potential and frictional 
profile for the above zero temperature case. The 
effective potential is scaled with respect to Vq. For 
the parameters chosen in the figure the current fiows 
in the negative direction as has been mentioned earlier. 



Following references [23, |2J| , one can obtain exact 
analytical expressions for the diffusion coefficient D as 



D 



j^:^'fu{x) 



(11) 



where I+{x) and I-{x) are as given below 
I+{x) ^ 



exp[i/;(a;)] / dy exp [- 7/;(y)](,12) 



I-{x) = exp[— 



x+L 



dy 



1 



A{y) 



exp [V^(y)](13) 



L here represents the period of the potential {— 2tt in 
our case). The Brownian particle takes a time r — 
L/v to traverse a distance L with a velocity v. The 
diffusive spread of the particle in the same time is given 
by < (Aq)2 >= 2Dt. The criterion to have a reliable 
transport is that the diffusive spread should be less 
compared to the distance traversed, i.e., < (Ag)^ 



2Dt < L"^. This in turn implies that Pe 
for coherent transport. 



Lv/D > 2 
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III. RESULTS AND DISCUSSIONS 

The velocity (v), diffusion constant (D) and the 
Peclet number (-Pe) are studied as a function of dif- 
ferent physical parameters. All the physical quantities 
are taken in dimensionless form. In particular, velocity 
and diffusion are normalized by (Vo/r?oi) and (Vb/ryo) 
respectively. Throughout our work we have set Vq 
and 7^0 to be unity. Similarly, T and T are scaled with 
respect to Vorjo and Vq respectively. We have used 
the globally adaptive scheme based on Gauss-Kronrod 
rules for numerical evaluations l25l. 
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FIG. 2: Plot of V, D and Pe vs (j) ior T = 0.36, A = 0.9 and 

r = o.o. 



In Fig. 12 we plot v, D and Pe function of 

phase difference at T = for a fixed noise strength 
of r = 0.36. The values of the physical parameters 
are mentioned in the figure caption. All the physical 
quantities are periodic function of 4> as expected and 
the velocity is zero at (j> = , tt and 27r. As can be 
seen from the plot of effective potential, Fig. the di- 
rection of current is negative for (j) upto tt and then it 
becomes positive. The current is antisymmetric around 
the point (p — tt, {V{(I)+tt) — —V{4>)). This is expected 
on general grounds for the case of a simple sinusoidal 
symmetric potential. The absolute value of current ex- 
hibits a maxima between to tt and tt to 2tt. The 
nature of currents being positive or negative can be 
readily inferred from the slope of the effective poten- 
tial (for example see Fig. . The diffusion coefficient 



is finite at all values of (f> and exhibits minima wher- 
ever the currents are zero {(j) — 0, tt, 2tt). However, the 
magnitude of minima is more at the point where cf) — tt. 
As expected, the diffusion coefficient is periodic in (p. 
However, it exhibits maxima at different values of 
than that for current. In the plot for Pe as a function 
of (j) wc have drawn a dotted line as a guide for eye. In 
regions where Pe > 2 the transport is said to be coher- 
ent. It is evident from the figure that there exists wide 
range of 4> in which the transport is coherent {Pe > 2) 
along with adjoining regions where the transport is less 
coherent {Pe < 2). These regions are however sensi- 
tively dependent on other physical parameters. In the 
present case of sinusoidal potential Peclct number as 
high as 4 can be obtained by properly tuning the pa- 
rameters. 




FIG. 3: Plot of V, D and Pe vs for F = 0.36, = 1.0, 
A = 0.9 and T = 0.0. 

To understand the role of spatial asymmetry in the 
potential we have included an asymmetry in the po- 
tential such that V{q) — Vo[sin(g) — ^/4sin(2g)] where 
fi is the asymmetry parameter. Fig. |21 shows the be- 
haviour of V, D and Pe for this simple asymmetric 
case at zero temperature with the other physical pa- 
rameters kept the same as in Fig. |2 The potential 
asymmetry parameter fi affects the behaviour of veloc- 
ity, diffusion and thereby Pe dramatically. We have 
considered the case of maximal asymmetry (/i — 1). 
We first notice that the simple symmetry observed for 
the case of periodic potential no more holds true in 
the presence of spatial asymmetry in potential. The 
magnitude of velocity is nonzero for values , tt and 
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27r. The velocity and diffusion constant in the region 
between to tt are suppressed while in the region be- 
tween TT and 27r are enhanced as compared with sym- 
metric case. Consequently, the presence of asymmetry 
can enhance or suppress the coherence in the transport. 
For the chosen parameters, Peclet number as high as 
5 is obtained. Moreover, for a given Pe either in- 
creases or decreases monotonically with the asymmetry 
parameter ^ (0 < /i < 1) which has been verified sep- 
arately. In our further analysis we restrict to the case 
of symmetric potential alone. 
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FIG. 4: Plot of 
and r = 0.01. 
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D and Pe vs F for 



0.627r, A = 0.9 



In Fig. 2] we plot velocity, diffusion and Pe as a func- 
tion of external noise strength F with cf) = 0.627r and 
A = 0.9 at finite temperature T — 0.01. It should be 
noted that velocity is negative in the entire range. The 
velocity is initially zero for F equal to zero and then 
increases with F and saturates to a constant value at 
higher values of noise strength. On the contrary, the 
diffusion constant keeps increasing monotonically with 
F. It is also clear that the external noise play a con- 
structive role in optimizing the coherence in transport 
i.e., the Pe exhibits a peak as a function of F. 

In Fig. [SI we plot velocity, diffusion and Pe as a 
function of temperature T with ~ 0.627r, A — 0.9 and 
F = 0.36. This value of noise strength corresponds to 
an optimal value for Pe in Fig. ^ The noise induced 
current is negative. The current exhibits a peak and 
then decreases to zero with temperature. This is 
expected because higher temperatures overshadows 
the effect of potential and frictional inhomogenieties 
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FIG. 5: Plot of 
and A = 0.9. 



Temperature 



D and Pe vs T for F = 0.36, (j> = 0.627r 




FIG. 6: Plot of V, D and Pe vs A for F = 0.36, = IA-k 
and r = 0.01. 



thereby suppressing the ratchet effect. In contrast, 
the diffusion constant increases monotonically. As 
opposed to the case of Fig. 0] we observe that internal 
fluctuations or temperature degrades the coherence 
in the transport i.e., Pe decreases with increase in 
temperature. 

In Fig. IHl we plot velocity, diffusion and Pe as 
a function of A, the amplitude of oscillation of the 
friction coefficient (0 < A < 1) for = 1.47r, F = 0.36 
and T = 0.01. All the physical quantities, namely, 
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velocity, diffusion and Pe increases monotonically 
with A. Thus the increase in A makes the transport 
more coherent. 



IV. CONCLUSIONS 

We have studied the coherence or reUabUty of trans- 
port of an overdamped Brownian particle in a frictional 
ratchet system with an underlying sinusoidal potential 
in the presence of external Gaussian white noise fluctu- 
ations. The frictional inhomogeneities along with ex- 
ternal fluctuations lead to a noise induced current or 
transport. The attained noise induced transport is al- 
ways accompanied by a diffusive spread which inturn 
makes the transport to be less reliable. We have shown 
that frictional inhomogeneities with respect to the un- 
derlying potential can make the transport coherent or 
incoherent. While the external noise (F) optimizes the 
coherence in transport the internal noise (T) degrades 
the coherence. 

In our present case, as mentioned in the beginning, 
the transport can be associated with an effective po- 
tential and an effective space dependent diffusion con- 



stant. The effective potential exhibits a tilt as a func- 
tion of system parameters (Fig.^. By looking at this 
effective potential one can infer only the direction of 
current and not its magnitude. The particle motion 
in this effective potential is determined by two time 
scales, (i) escape from the potential minima over the 
barrier along the effective bias followed by (ii) the re- 
laxation into next minima. The coherent transport in 
homogeneous medium is obtained when the relaxation 
time dominates the transport 0, 0| as conipared with 
the escape time. For details see reference We 
would like to emphasize that the dynamics of the par- 
ticle in our present problem arises due to the complex 
interplay between the potential, internal fluctuations, 
frictional proflle and external fluctuations. This is am- 
ply reflected in the fact that both the effective potential 
and space dependent diffusion constant change their 
nature as we change the system parameters. Thus apri- 
ori analysis or prediction of the behaviour of the sys- 
tem in regard to transport coherence is a difficult task. 
It may not be surprising that by choosing appropriate 
asymmetric potential and frictional profile along with 
other parameters one may obtain much higher trans- 
port coherence as observed in 0. 
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